Tuning surface metallicity and ferromagnetism by hydrogen adsorption at the polar 

ZnO(OOOl) surface 
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The adsorption of hydrogen on the polar Zn-ended ZnO(OOOl) surface has been investigated by 
density functional ab-initio calculations. An on top H(l x 1) ordered overlayer with genuine H-Zn 
chemical bonds is shown to be energetically favorable. The H covered surface is metallic and spin- 
polarized, with a noticeable magnetic moment at the surface region. Lower hydrogen coverages lead 
to strengthening of the H-Zn bonds, corrugation of the surface layer and to an insulating surface. 
Our results explain experimental observations of hydrogen adsorption on this surface, and not only 
predict a metal-insulator transition, but primarily provide a method to reversible switch surface 
magnetism by varying the hydrogen density on the surface. 

PACS numbers: 71.20.Tx, 73.22.-f, 71.70.Ej 



ZnO is one of the most technologically important metal 
oxides, which holds great promise for applications as 
a wide band gap semiconductor: it shows the largest 
charge-carrier mobility among oxides, extraordinary cat- 
alytic properties and the unusual coexistence of trans- 
parency and conductivity The recently reported 
high temperature (HT) ferromagnetism in ZnO-based 
thin layers and nanostructures [2J turns it additionally 
into a potential HT semiconducting ferromagnet, which 
could be used in magnetoelectric and magnetotransport 
devices tuning simultaneously charge and spin How- 
ever, at present, and despite numerous experimental and 
theoretical studies, the mechanism behind the HT mag- 
netic order in ZnO is still under debate 0, H, @, 0] • 

Particular attention deserves the role of H doping in 
ZnO. Hydrogen is a very reactive element that exhibits 
qualitatively different behaviour in different media. It 
is amphoteric, can act either as donor (H+) or acceptor 
(H-), can occupy different lattice sites and even modify 
the host structure, and in general counteracts the con- 
ductivity of the host. Isolated hydrogen has been found 
to act as a shallow donor in bulk ZnO, and thus it has 
been attributed as a source of the unintentional n-type 
conductivity exhibited by ZnO [8j. 

At surfaces, the interaction with ambient H is al- 
most unavoidable even in ultra-high- vacuum conditions. 
Moreover, the presence of hydrogen has a pronounced 
influence which can even change the surface electronic 
properties. Understanding the interaction of hydrogen 
with the ZnO surfaces is crucial in order to control 
the properties of ZnO-based low-dimensional structures. 
In recent years, systematic investigations of the differ- 
ent ZnO low- indexed surfaces have been performed Q. 
There is agreement in the formation of an ordered hy- 
drogen overlayer on both the non-polar ZnO (1010) and 
the polar O-ended ZnO (0001) surfaces. Much few work 
has been devoted to the Zn-terminated ZnO (0001) sur- 
face, where the interaction of H atoms is thought to be 



the weakest among the ZnO surfaces, since the bind- 
ing energy of ZnH pairs should be considerably weaker 
than that of OH bonds. An experimental study revealed 
that, exposing this surface to atomic hydrogen, an or- 
dered (1 x _1) overlayer consisting of Zn- hydride species 
is formed [lfj. However, larger exposures to both atomic 
and molecular hydrogen cause the unexpected unstabil- 
ity of the H (1 x 1) overlayer, leading to a complete loss 
of lateral order which reflects a random distribution of H 
adatoms. 

The aim of this Letter is to show that the adsorption of 
atomic hydrogen on the polar Zn-ZnO(OOOl) surface gives 
rise to an ontop H(l x 1) ordered overlayer with genuine 
H-Zn chemical bonds. As shown below, surface mag- 
netism and metallicity are distinct characteristics of this 
two-dimensional (2D) H(l x 1) ordered overlayer. Fur- 
thermore, the partial coverage of the surface with H leads 
to reinforcement of the H-Zn bonds, corrugation of the 
surface layer and, more interestingly, the emergence of 
a spin-paired insulating state. Hence, a metal-insulator 
transition accompanied by the extinction of the magne- 
tization can be driven by reducing the H density on the 
surface from 1 ML down to 1/2 ML. Further decrease of 
the H coverage restores the surface metallicity. 

Our calculations are based on density functional the- 
ory employing norm-conserving pseudopotentials and lo- 
calized numerical atomic orbitals (AO) as implemented 
in the SIESTA code[llj. We consider both the local spin 
density approximation (LSDA) and the generalized gra- 
dient approximation (GGA), with the same basis set and 
parameters of Ref . 5] . Hydrogen is described by a Dou- 
ble Z Is AO. In order to ensure that our conclusions 
are robust against the choice of the exchange-correlation 
functional, we have performed additional total energy op- 
timizations with both the PBE0 and HSE hybrid func- 
tional [12, EH using the Projector Augmented Wave 
(PAW) method implemented in the VASP code [3]. We 
employed similar parameters to those used in Ref. [HI], 
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where such functionals were sucessfully applied to study 
the defect energetics in ZnO. 

The surfaces are modeled by periodically repeated 
slabs containing between 15 and 17 atomic planes, sep- 
arated by a vacuum region of at least 20 A. Bulk-like 
behaviour is always attained at the innermost central lay- 
ers. We calculate both symmetric and asymmetric slabs 
about the central plane, though we are always subject to 
the lack of inversion symmetry of the wurtzite structure. 
We use (1 x 1) and (2 x 2) 2D-unit cells to model the sur- 
faces with adsorbed hydrogen. All the atomic positions 
are allowed to relax until the forces on the atoms are less 
than 0.03 eV/A. Brillouin Zone integrations have been 
performed on a 12 x 12 x 1 Monkhorst-Pack supercell 
(12 x 12 x 8 for bulk structures). Careful convergence in 
the k-mesh and the inclusion of relaxations are essential 
in order to accurately describe the adsorption states and 
magnetism. 
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FIG. 1: PDOS projected on the H (filled curve) and neigh- 
bouring Zn and O atoms of the (top) bare Zn-Zn(OOOl) sur- 
face, (middle) 1 ML H covered surface, and (bottom) bulk 
ZnO with a H impurity in the MBC. Left panels correspond 
to LDA, right ones to HSE. 

For 1 ML coverage, the stable adsorption site of H 
was determined to be atop the Zn atoms after exhaustive 
minimization considering several surface and subsurface 
positions, including fee and hep hollow, bridge and off- 
symmetric sites. We present in Figure Q] the density of 
states projected (PDOS) onto the outermost layers for 
the bare and H covered surface, both using LDA and 
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FIG. 2: Side view of the spin density distribution at the 
ZnO(OOOl) surface covered by 1 ML H. 



the HSE functional. For comparison, the PDOS corre- 
sponding to a substitutional H atom at an oxygen site in 
bulk ZnO -multicenter bonds configuration (MBC) [8J- is 
also shown. As expected, the HSE functional provides a 
better description of the band gap and increases the lo- 
calization of the Zn d-states. The calculated gaps using 
LDA and HSE for bulk ZnO are 0.80 and 3.45 eV, respec- 
tively, compared to the experimental value of 3.60 eV. 
Nevertheless, both calculations provide the same phys- 
ical picture. The clean surface exhibits the well-known 
metallic character, with surface states (SS) at the bottom 
of the conduction band (CB) [nij]. However, H adsorp- 
tion creates hybridized bonding orbitals at the valence 
band (VB) edge, depleting the Zn 4s SS at the CB. The 
Fermi level lies in the H derived states and thus the sur- 
face is metallic and p-doped. This is opposite to the bulk 
MBC, for which the H-bonding state lies deep below the 
VB with a total charge not large enough to completely 
deplete the Zn derived CB states. 

A distinct characteristic of the H covered surface in the 
atop geometry is the spin-polarization of the bands. It is 
not restricted to the H layer, but extends into the ZnO 
subsurface leading to a net magnetization of the surface 
region. This is clearly seen in Fig. [21 The magnetic 
moments at the H, Zn and O layers are 0.29, 0.09 and 
0.09 /is, respectively, and the magnetic energy gain is 
around 70 meV. Similar values are obtained in the HSE 
calculation, the total magnetization differing in less than 
10%. Summarizing, H adsorption at the ML coverage 
gives rise to a p-doped metallic surface with the Fermi 
level pinned at the H-derived bands and a net surface 
magnetic moment of around 0.5 /ib- 

In order to understand the emergence of the spin po- 
larization when going from the bare Zn surface to the 
complete H overlayer, we have modelled partial Hydrogen 
coverages of 3/4, 1/2 and 1/4 MLs with LDA. The cor- 
responding geometric and electronic structures are dis- 
played in Figure[3]and further details are given in Tabled] 
The adsorption of H creates low dispersion states at the 
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FIG. 3: (Top) Side views of the ZnO(OOOl) surface for decreasing H coverage from left to right. Oyxgen, Zn and H atoms are 
colored in red, blue and white, respectively. Below each sketch we provide for each case the corresponding: (middle) PDOS 
projected on the first surface layer and resolved in fc-space[T3|, the BZ for the (lxl) and (2x2) 2D cells being depicted at the 
right of the bottom panel; (bottom) H-Zn and H-0 (magnified by a factor of 5) COOPs. 



top of the VB, which are spin-splitted for 1 ML coverage 
(with only one spin component completely filled), and 
start to merge for lower coverages until they become de- 
generate for 1/2 ML of H. The main result is that this 
loss of spin-polarization is accompanied by a metal to 
insulator transition of the surface. Thus, while at 3/4 
ML the surface is still ferromagnetic and metallic, for 
a coverage of 1/2 ML the ground state corresponds to a 
non-magnetic insulator with the Fermi level above the hy- 
drogen derived band. Further decrease of the hydrogen 
coverage results in partial occupation of the CB states 
due to the increased number of unsaturated Zn dangling 
bonds, restoring the surface metallicity. Therefore, a 
metal to insulator transition can be tuned by varying the 
H-coverage back and forth. Even more, by varying the 
Hydrogen density reversible switch of surface magnetism 
can be achieved. 

A deeper insight about this remarkable phenomenon 
can be obtained regarding the nature of the bonds -see 
Table[l] In all cases, the H-Zn bond length is significantly 
shorter than at the bulk (~ 2 A in the MBC), a hint of the 
formation of genuine and stronger chemical bonds at the 
surface. Furthermore, the bond length reduces as the H 
coverage diminishes, evidencing a reinforcement of the H- 
Zn bonds. This bond strengthening is consistent with the 
variation of the Mulliken charges and the corresponding 



TABLE I: H-related surface energy (E stlr /, in eV), bond dis- 
tances (&ab, in A), and Mulliken charges (Qa) of the surface 
atoms for different H coverages and for the bare Zn-ZnO(OOOl) 
surface. The last column refers to undoped bulk ZnO. 

1H 3/4 H 1/2 H 1/4 H bare bulk" 

d H -zn 1.62 1.60 1.57 1.56 

dzn-p 0.71 0.78/0.21 0.91/0.20 1.02/0.14 0.39 0.61 
Qh 1.09 1.12 1.17 1.21 

Qzn 11-36 11.38/11.25 11.42/11.24 11.45/11.27 11.37 11.23 
E aurf -1.53 -1.86 -2.51 -3.22 ~ 



crystal overlap populations (COOPs) at the lower panel 
of Figure [31 The charge at both the H and the Zn in- 
creases as the H coverage reduces, indicating that a larger 
amount of charge is shared between the two atoms. The 
COOPs confirm this scenario, their positive value corre- 
sponding to bonding states [l8(, which are progressively 
filled as the H coverage reduces. Noticeably, Hydrogen 
also interacts with the O in the subsurface layer forming 
an antibonding state. 

The reinforcement of the H-Zn bonds can be under- 
stood regarding the atomic structures in Figure [3] and 
Table [fl For the clean surface we find a contraction 
of the first double layer spacing, dz n -o, m agreement 
with previous calculations [l6|. Contrary, the surface 
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completely covered with H exhibits a slightly expanded 
Zn-0 distance. For partial H coverages the surface Zn 
atoms become inequivalent, leading to two different Zn- 
O interlayer distances and a large corrugation of the Zn 
plane: those bonded to H experience an outward relax- 
ation while the unbonded ones show an inward relaxation 
even stronger than that at the bare surface, becoming al- 
most coplanar with the oxygen plane. For deeper layers 
little relaxation is found. There is a strong correlation 
between the Zn-0 interlayer distance and the redistri- 
bution of charge: Zn atoms bonded to H -large dz n -0- 
show an increase in their Mulliken charges, while those 
not bonded -small dz n -o- lose charge approaching the 
bulk values. 

The energies provided in Table [J for the H covered 
surfaces also support the formation of strong bonds with 
a significant covalcnt character. They are calculated as 
the difference between the total energy of the relaxed slab 
and those of the relaxed isolated ZnO and hydrogen slabs. 
Their negative values indicate that the chemisorption of 
hydrogen on the Zn(0001) surface is exothermic. In fact, 
at a hydrogen density of 1/2 ML the surface energy ex- 
ceeds by 140 eV the formation of the H2 molecule in the 
gas phase, 2.37 eV/H, making adsorbed hydrogen highly 
stable against desorption. In addition, there is a net re- 
duction of the surface energy for all H coverages with 
respect to the bare surface. This is specially significant 
for partial H coverages, where the existence of inequiva- 
lent Zn sites allows for the corrugation of the Zn layer and 
the subsequent decrease of surface energy. These results 
consistently explain experimental observations, not only 
about the formation of a H(l x 1) ordered overlayer, but 
also the loss of surface order under prolongated hydrogen 
exposures fioj ]. As shown above, partial H adsorption is 
energetically more favorable than the bare surface, which 
could explain the experimentally observed enhancement 
of the chemical reactivity of surfaces previously exposed 
to H with respect to the pristine ones. Also the almost 
coplanar positions of the Zn and O atoms in the bare 
regions of the partially H covered surface can account for 
the oxygen signature observed in the XPS spectra [icj |. 

An additional important conclusion can be extracted 
with respect to the ferromagnetism measured in ZnO, 
both undoped and doped with magnetic impurities. It 
has been observed to vary with the oxygen partial pres- 
sure, and consequently its existence has been associated 
to oxygen vacancies. However, calculations as well as 
careful experiments indicate that oxygen vacancies can- 
not be a source of magnetism, but instead unsaturated 
oxygen either due to Zn vacancies or to O-terminated 
surfaces 0, d, 0, [HJ. Our results strongly suggest that 
unintentionally adsorbed hydrogen may take active part 
in the observed magnetism, particularly for undoped ZnO 
nanocrystals. Adsorbed hydrogen is sensitive to the oxy- 
gen partial pressure, since the probability of formation 
of HO complexes increases with the oxygen chemical po- 



tential, leading to the subsequent extinction of the mag- 
netism reported here. Moreover, the hydrogen electron 
spin has been measured in In 2 03 whe n hy drogen acts as 
an oxygen vacancy passivating center [20(. 

As a final remark, the Zn-ZnO(OOOl) surface may 
undertake a triangularly shaped reconstruction under 



specific growth conditions 21]. We have investigated 



the effect of such reconstruction through calculations of 
stepped surfaces, and the general conclusions are anal- 
ogous to those presented here for the unreconstructed 
surface, although strong O-H bonds are also formed at 
step edges. 

In summary, atomic hydrogen adsorbs on the Zn- 
ZnO(OOOl) polar surface atop the Zn atoms, forming 
strong H-Zn bonds and leading to a metallic, p-doped 
surface with a net magnetic moment. As the H cover- 
age diminishes, there is a reinforcement of the remaining 
H-Zn bonds accompanied by the removal of the valence 
band holes and the subsequent extinction of magnetism. 
We predict that controlling the Hydrogen coverage can 
serve to tune a metal-insulator transition and to achieve 
reversible switch of surface magnetism. 
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